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Abstract. We study discrete complex analysis and potential theory on a large fam- 
ily of planar graphs, the so-called isoradial ones. Along with discrete analogues of 
several classical results, we prove uniform convergence of discrete harmonic measures, 
Green's functions and Poisson kernels to their continuous counterparts. Among other 
applications, the results can be used to establish universality of the critical Ising and 
other lattice models. 



1. Introduction 

1.1. Motivation. This paper is concerned with discrete versions of complex analysis 
and potential theory in the complex plane. There are many discretizations of harmonic 
and holomorphic functions, which have a long history. Besides proving discrete ana- 
logues of the usual complex analysis theorems, one can ask to which extent discrete 
objects approximate their continuous counterparts. This can be used to give "discrete" 
proofs of continuous theorems (see, e.g., |L-F55] for such a proof of the Riemann map- 
ping theorem) or to prove convergence of discrete objects to continuous ones. One of 
the goals of our paper is to provide tools for establishing convergence of critical 2D 
lattice models to conformally invariant scaling limits. 

There are no "canonical" discretizations of Laplace and Cauchy-Riemann opera- 
tors, the most studied ones (and perhaps the most convenient) are for the square grid. 
There are also definitions for other regular lattices, as well as generalizations to larger 
families of embedded into C planar graphs. 

We will work with isoradial graphs (or, equivalently, rhombic lattices) where all 
faces can be inscribed into circles of equal radii. Rhombic lattices were introduced by 
R. J. Duffin [Duf68j in late sixties as (perhaps) the largest family of graphs for which 
the Cauchy-Riemann operator admits a nice discretization, similar to that for the 
square lattice. They reappeared recently as isoradial graphs in the work of Ch. Mercat 
[MerOlj and R. Kenyon [Ken02j . as the largest family of graphs were certain 2D statisti- 
cal mechanical models (notably the Ising and dimer models) preserve some integrability 
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properties. Note that isoradial graphs can be quite irregular - see e.g. Fig. [T]A.. It was 
shown by R. Kenyon and J.-M. Schlenker |KSch04j that many planar graphs admit 
isoradial embeddings - in fact, there are only two topological obstructions. Also isora- 
dial graphs have a well-defined mesh size 8 - the common radius of the circumscribed 
circles. 

It is thus natural to consider this family of graphs in the context of universality 
for 2D models with (conjecturally) conformally invariant scaling limits (as the mesh 
tends to zero). 

The primary goal of our paper is to provide a "toolbox" of discrete versions of 
continuous results (particularly "hard" estimates) sufficient to perform a passage to 
the scaling limit. Of particular interest to us is the critical Ising model, and this 
paper starts a series devoted to its universality (which means that the scaling limit is 
independent of the shape of the lattice). We plan to establish convergence of certain 
discrete holomorphic observables, and then deduce convergence of interfaces to Oded 
Schramm's SLE curves and corresponding loop ensembles. See [ChSm08j, [Sm06j for 
the strategy of our proof and [Sm07| for the square lattice case. 

Our results can also be applied to other lattice models. The uniform convergence 
of the discrete Poisson kernel (11.31) already implies universality for the loop-erased ran- 
dom walks on isoradial graphs. Namely, our paper together with |LSchW04] implies 
that their trajectories converge to SLE(2) curves (see Sect. 3.2, especially Remark 3.6, 
in |LSchW04] ). There are several other fields where discrete harmonic and discrete 
holomorphic functions defined on isoradial graphs play essential role and hence where 
our results may be useful: approximation of conformal maps |Biick08] ; discrete inte- 
grable systems [BMS05J; and the theory of discrete Riemann surfaces |Mer07] . 

Local convergence of discrete harmonic (holomorphic) functions to continuous 
harmonic (holomorphic) functions is a rather simple fact. Moreover, it was shown by 
Ch. Mercat [Mer02j that each continuous holomorphic function can be approximated 
by discrete ones. Thus, the discrete theory is close to the continuous theory "locally." 
Nevertheless, until recently almost nothing was known about the "global" convergence 
of the functions defined in discrete domains as the solutions of some discrete boundary 
value problems to their continuous counterparts. This setup goes back to the seminal 
paper by R. Courant, K. Friedrichs and H. Lewy [CFL28j . where convergence is estab- 
lished for harmonic functions with smooth Dirichlet boundary conditions in smooth 
domains, discretized by the square lattice, but not much progress has occurred since. 
For us it is important to consider discrete domains with possibly very rough boundaries 
and to establish convergence without any regularity assumptions about them. Besides 
being of independent interest, this is indispensable for establishing convergence to Oded 
Schramm's SLEs, since the latter curves are fractal. 

1.2. Preliminary definitions. The planar graph T embedded in C is called isoradial 
iff each face is inscribed into a circle of a common radius 5. If all circle centers are 
inside the corresponding faces, then one can naturally embed the dual graph T* in 
C isoradially with the same 5, taking the circle centers as vertices of T*. The name 
rhombic lattice is due to the fact that all quadrilateral faces of the corresponding 
bipartite graph A (having the vertex set TUT*) are rhombi with sides of length S (see 
Fig.fTK). In our paper we also require the following mild but indispensable assumption: 
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(A) (B) 

Figure 1. (A) An isoradial graph T (black vertices, solid lines), its 
dual isoradial graph T* (gray vertices, dashed lines), the corresponding 
rhombic lattice or quad-graph (vertices A = T U T*, thin lines) and the 
set <0 = A* (white diamond-shaped vertices). (B) Local notations near 

u e r. 



(4k) the rhombi angles are uniformly bounded away from and it (in 
other words, all these angles belong to [rj, tt — i]] for some fixed 77 > 0). 

We will often work with rhombi half-angles, denoted by 6, for which we then have 

7T 71 

< const < 9 < const < — . 

2 2 

Note that condition (4&) implies that for each u\, 112 € V the Euclidean distance \u2— u\\ 
and the combinatorial distance 5 ■ dr(ui,U2) (where dr(ui,U2) is the minimal number 
of vertices in the path connecting u\ and u<i in T) are comparable. Below we often use 
the notation const for absolute positive constants that does not depend on the mesh 5 
or the graph structure but, in principle, may depend on 77. 

The function H : fif — >• R defined on some subset (discrete domain) f2p of T is 
called discrete harmonic, if 

n 

^2time s -(H(u s )-H(u)) = (1.1) 

s=l 

at all u e fip where the left-hand side makes sense. Here 9 S denotes the half-angles 
of the corresponding rhombi, see also Fig. QJ3 for notations. As usual, this definition 
is closely related to the random walk on V such that the probability to make the next 
step from u to Uk is proportional to tan^.. Namely, RW(t+l) = RW(t) +^Rwctl' where 
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the increments £w are independent with distributions 

P(£ M = u k -u) = „„ an 6lfc for A; = 1, ..,n. 
2^ s= i tan S 

Under our assumption all these probabilities are uniformly bounded from 0. Note that 
the choice of tan# s as the edge weights in (11. ip gives 

*rp c 1 W |t c 1 n , E [( Re ^) 2 ] = E[(Im^) 2 ] = T w 
E[Re£ u ] = E[Im&J =0 and (1.2) 

E[Ref u Imf u ] = 0, 

where T u = <5 2 • X^=i sin 29 S /Yl™=i tan S (see Lemma I2T2}) . Our results may be directly 
interpreted as the convergence of the hitting probabilities for this random walk. More- 
over, condition (4) implies that quadratic variations satisfy < const 5 2 ^ T u ^ 5 2 , 
so after a time re-parameterization according to (11.21) . this random walk converges to 
standard 2D Brownian motion. 

1.3. Main results. Let Of C T be some bounded, simply connected discrete domain 
and Int Of, <90f denote the sets of interior and boundary vertices, respectively (see 
Sect. 12.11 for more accurate definitions). For u G Int Of and E C <90f the discrete 
harmonic measure uj 5 (u; E; Of) is the probability of the event that the random walk 
on T starting at u first exits Of through E. Equivalently, o/(-; E; Of ) is the unique 
discrete harmonic function in Of having boundary values 1 on E and on <90f \ E. 

We prove uniform (with respect to the shape Of and the structure of the under- 
lying isoradial graph) convergence of the basic objects of the discrete potential theory 
to their continuous counterparts. Namely, we consider 

• discrete harmonic measure u/(-; a s b s ; Of) of arcs a s b s C <90f (Sect. 13.31) ; 

• the discrete Green's function G s n6 (-;v s ), v s G Int Of (Sect. 13.31) ; 

• the discrete Poisson kernel 

normalized at the interior point v 5 G Int Of (Sect. 13.41) ; 

• discrete Poisson kernel P s s { 4 ] o 5 ; Of), a s G 90f , normalized at the boundary by 
the discrete analogue of the condition -§-P\ s = —1 (we assume that Of has the 
"straight" boundary near o s G <90f , see the precise definitions in Sect. 13.51) . 

Moreover, we also prove uniform convergence for the discrete gradients of these func- 
tions (which are discrete holomorphic functions defined on subsets of <^ = A*, see Sect. 
12.41 and Definition 13.61 for further details). 

1.4. Organization of the paper. We begin with the exposition of basic facts con- 
cerning discrete harmonic and discrete holomorphic functions on isoradial graphs. The 
larger part of Sect. M follows |Duf68j . jMerOl] . [Ken02j . |Mer07j and |Buck08j . Unfor- 
tunately, none of these papers contains all the preliminaries that we need. Besides, 
the basic notation (sign and normalization of the Laplacian, definition of the discrete 
exponentials and so on) varies from source to source, so for the convenience of the 
reader we collected all preliminaries in the same place. Note that our notation (e.g., 
the additive normalization of the discrete Green's functions) is chosen to be as close in 
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the limit to the standard continuous objects as possible. Also, we prefer to deal with 
functions rather than to use the language of forms or cochains |Mer07] which is more 
adapted for the topologically nontrivial cases. 

The main part of our paper is Sect. [3J where the convergence theorems are proved. 
The proofs essentially use compactness arguments, so it does not give any estimate for 
the convergence rate. Thus, as in |Sm07j , we derive the "uniform" convergence from the 
"pointwise" one, using the compactness of the set of bounded simply connected domains 
in the Carathedory topology (see Proposition 13.71) . The other ingredients are the 
classical Arzela-Ascoli theorem, which allows us to choose a convergent subsequence of 
discrete harmonic functions (see Proposition EH]) and the weak Beurling-type estimate 
(Proposition 12.101) which we use in order to identify the boundary values of the limiting 
harmonic function. We prove C7 1 -convergence, but stop short of discussing the C°° 
topology since there is no straightforward definition of the second discrete derivative 
for functions on isoradial graphs (see Sect. 12.51) . Note however that a way to overcome 
this difficulty was sugg ested in [Biick08j . 

Acknowledgments. We would like to thank Vincent Beffara for many helpful com- 
ments. Parts of this paper were written at the Mathematisches Forschungsinstitut 
Oberwolfach during the Oberwolfach-Leibniz Fellowship of the first author. The au- 
thors are grateful to the MFO for the hospitality. 

2. Discrete harmonic and holomorphic functions. Basic facts 

2.1. Basic definitions. Approximation property. Let r = T 5 be some (infinite) 
isoradial graph embedded into C. Let fl s = Fqs U Eqs U Vqs, where Fqs, Eqs and 
Vqs are sets of (open) faces, (open) edges and vertices of T 5 (see Fig. [2]A.). If Q s C C 
is a (simply) connected domain (i.e., open subset of C) and the skeleton E n s U Vqs 
is connected, we call fl 5 the continuous counterpart of the (simply) connected 
discrete domain f2 r . The latter is defined as the set of vertices f2f = Int f2f U dfl^, 
where the sets of interior and boundary vertices are given by 

Int tt 5 r := V Q s and dtt s r := {(a; a int ) : (a int a) G E n s, a int G Vqs, a V n s}, 

respectively. The reason for a more complicated definition of a boundary vertex (as a 
pair of an interior and an exterior vertices) is that it corresponds to a prime end in the 
continuous setting. In particular, the same vertex will play the role of several boundary 
vertices, when it can be approached from different directions inside the domain Q 5 - 
see e.g. vertices b and c in the Fig. [2JA). However, when no confusion arises we will 
often denote boundary vertices by a, not indicating explicitly the edge (a int a) that ends 
at a. 

For a domain HcCwe define its discretization il 5 (and, further, f2 r ) by taking 
the (largest) connected component of T lying inside fl as Vqs, and the set of all edges 
and faces incident to Vqs as Eqs and Fqs, respectively (see Fig. EB, Fig. E]A.). Note 
that, by definition, <90f n Q = 0. 

Further, we introduce the weight of a vertex u G V by 




(2.1) 
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Figure 2. (A) Discrete domain. The interior vertices are gray, the 
boundary vertices are black and the outer vertices are white. Both b and c 
have two interior neighbors, and so we treat ftW = (6; b^ t ), = (6; b^ t ) 
and c^ l \ as different elements of <9fi r . (B) Discrete half-plane HP 5 
and discrete rectangle R S (S,T). The lower, upper and vertical parts of 
dR^{S,T) are denoted by L s r (S), U^{S,T) and V*(S, T), respectively. 



where 9 S are the half-angles of the corresponding rhombi (see Fig. H]B). An easy exercise 
in trigonometry shows that /if(w) * s the area of the dual face W(u) = Wiw 2 ..w n , 
u ~ w s g r*. 

Let : Q 5 — > C be a Lipschitz (i.e., satisfying |0(wi) — 0(«2)| ^ — w 2 |) 

function and <p 5 = </>| n i be its restriction to fi r . Note that all points in the dual face 
W(u) are 5-close to u, so approximating values of (f> on W(iz) by its value at u we arrive 
at the inequality 



\J 4> s (u)^(u) + (j) S {a)jJ n s (a) — // 4>{x+iy)dxdy 



< C5 • Area(0 5 ) (2.2) 



with the same constant C . Here /x^ 4 (a) is the area of W(a) n fi* 5 (more accurately, the 
area of the connected component of W (a) D intersecting with the edge (aai nt )). 

Definition 2.1. Lei fif 5e some connected discrete domain and H : f2 r — > R. VFe 
define the discrete Laplacian of H at u G Int f2 r as 

[A^u) := * Y, ^6 S ■ [H{u s )-H{u)). 
VFe caZ/ H discrete harmonic in VL 5 T iff [/S. 5 H] {u) = at all interior points u G Int fi r . 
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It is easy to see that discrete harmonic functions satisfy the maximum principle: 

max H(u) = max H(a). 

Furthermore, the discrete Green's formula 

[GA s H-HA s G](u)4(u) = ^ tznd aaint ■ [H(a)G(a int ) - H(a int )G(a)] 

holds true for any two functions G, H : Of — > R. Here aa . nt denotes the half-angle of 
the rhombus having; diagonal. 

Lemma 2.2 (approximation property). Let G C 3 be a smooth function defined 

in the disc B(u, 25) C C for some u G T. Denote by <p s its restriction to T. Then 

(i) A s <p s = for all linear and A s <p s = A<p = 2(a + c) ; if <fi(x + iy) = ax 2 + bxy + cy 2 . 

(ii) 

I [A^K") - [A0](«) I < const -5 • max|L) 3 0|. 

Proof. We start by enumerating neighbors of « as Mi,...,«„ and its neighbors on the 
dual lattice as wi, . . . , w n - see Fig. [IB)- Obviously, A 5 <p 5 = 0, if <fi is a constant. Since 

^ tan6» s ■ (u a —u) = -i ^ (w s+ i-u7 s ) = 0, 

one obtains A* 5 */)* 5 = for linear functions x = Rew and y = Ima. Similarly, 

tan6» s ■ (u 2 s -u 2 ) = -i ^2(w s+1 -w s )(u+u s ) = -i ^2(w 2 +1 -w 2 ) = 0, 

U S ^U U S ^U U S ^U 

so A S (f) s = for x 2 —y 2 = Keu 2 and 2xy = Imu 2 . The result for x 2 +y 2 follows from 

tan °s ■ \u s -u\ 2 = 25 2 Y sin2^ s = 4//f(u), 

thus proving (i). Finally, the Taylor formula implies (ii). □ 
2.2. Green's function. Dirichlet problem. Harnack lemma. Lipschitzness. 

Definition 2.3. Let u G T. We call H = C7 r (-; u ) : T — > R £/ie free Green function 

z/f ^ satisfies the following: 

(i) [A 5 H}(u) = for all u ^ u and ii s r (u ) ■ [A 5 H](u ) = 1; 

H(u) = o(\u — uq\) as \u — Uq\ — > oo; 
(mj H{uq) = 2^(log5— 7Euier — log2) ; where 7E U ier i/ie Euler constant. 

Remark 2.4. VFe use £/ie nonstandard normalization at uq (usually the additive con- 
stant is chosen such that G(uq;uq) = 0) in order to have convergence to the standard 
continuous Green's function ^log \u — Uq\ as the mesh 5 goes to zero. 

Theorem 2.5 (improved Kenyon's theorem). There exists a unique Green's func- 
tion Gt(-;uo). Moreover, it satisfies 

Gr(u] Uq) = — — log I it ito| + O ( - — - — -], M^tt , (2.3) 
2it \\ u ~ u o\ J 

uniformly with respect to the shape of the isoradial graph F and uq G T. 
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Proof. This asymptotic form for isoradial graphs was first obtained in |Ken02j . Some 
small improvements (the correct additive constant and the order of the remainder) were 
done in [Biick08j under the so-called quasicrystallic assumption (when the number of 
edge slopes in the corresponding rhombic embedding is finite). Actually, Kenyon's 
beautiful proof works well without any appeal to this assumption. We give the sketch 
of the proof (together with a slight but necessary modification) in Appendix IA.1I □ 

Let Q s r be some bounded connected discrete domain. It is well known that for 
each / : <9fif — > R there exists a unique discrete harmonic function H in fif such 
that H\qq^ = f (e.g., H minimizes the corresponding Dirichlet energy, see [Duf68j). 
Clearly, H depends on / linearly, so 

H(u)= u S (u;{a};n 5 T )-f(a) 

for all u G f2 r , where u s {u; •; f2 r ) is some probabilistic measure on <9f2p which is called 
harmonic measure at u. It is harmonic as a function of u and has a standard 
interpretation as the exit probability for the random walk on V (i.e. the measure of a 
set A C <9fif is the probability that the random walk started from u exits fif through 
A). 

Definition 2.6. Let u G Int f2f . We call H = G n s (•; uq) the Green function in f2f 

iff 

(i) [A s H](u) = for all u £ Int fif \ {u } and //f (u ) ■ \A 5 H}(u ) = I; 

(ii) H = on the boundary dQ^. 

Note that these two properties determine G n s(-;uo) uniquely. Namely, G aS = Gr — 
G* nS , where 

G* n ^ u o)= Yl ^(•;{4;^r)-G ! r(a;«o) (2.4) 

is the unique solution of the discrete boundary value problem 

A s H* = in Q 5 V , H* = G(-;u ) on dtt 5 r . 
Applying Green's formula to oj 5 {-; {a}; fif) and G n ^ (•; u ), one obtains 

lo\uq\ {a};Q 5 T ) = -tan6» aQ . nt • G n * (a int ; u ), where a = (a; a int ) G dQ^. (2.5) 

It was noted by U. Bucking |Biick08] that, since the remainder in (12.31) is of order 
0(5 2 \u — uq\~ 2 ), one can use R. Duffm's ideas [Duf53j in order to derive the Harnack 
Lemma for discrete harmonic functions. Recall that we denote by B^(uq,R) C T the 
discretization of the disc B(u ,R) C C. 

Proposition 2.7 (Harnack Lemma). Let u G T and H : B^(u , R) — > R be a 
nonnegative discrete harmonic function, 
(i) If Ui ~ Uq, then 

\H(u x ) - H(u )\ «C const . 6H( ^°\ 
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(ii) Ifui,u 2 G B^(u ,r) C Int B^(uq, R), then 

r 

const • — . 

R — r 

Proof. In order to make our presentation complete, we recall briefly the arguments 
from [Buck08j and [Duf53] in Appendix IA.21 □ 

Corollary 2.8 (Lipschitzness of discrete harmonic functions). Let H be discrete 
harmonic in Bp(u ,R) andui,u 2 G Bf,(u ,r) C Int B^(u , R) . Then 

\H{u2) — H(ui)\ ^ const — — , where M— max \H(u)\. 

R — r B 5 r (u ,R) 

Proof. By assumption (4) we can find a path u\ = VQV\V 2 ---Vk-iVk = u 2 , connecting 
u\ and u 2 inside B^(uo,r), such that k ^ const •|m2 — mi|5 _1 . Since ^ H + M ^ 2M, 
applying Harnack's inequality to H + M, one gets 

k— 1 

\H(u 2 ) - H{u x )\ ^ V \H(v j+1 ) - H( Vj )\ < const .\^A . □ 
z — ' o K — r 

3=0 

2.3. Weak Beurling-type estimates. The following simple fact is based on the ap- 
proximation property (Lemma 12.21) for the discrete Laplacian on isoradial graphs. 

Lemma 2.9. Let u G T, r>0 and B^(u ,r) be the discretization of the disc B(u ,r). 
Let a, b G dB^(uo, r) be two boundary points such that 

arg(6-w ) - arg(a-Mo) ^ 

Then, 

u 5 (u; ab; B[(0, r)) ^ const > for all u G (u , |r), 
where ab denotes the discrete counter clockwise arc from a to b (see Fig. 

Proof. Fix some smooth function O : -8(0, |) — > R such that 

(ia) 0o(C) ^ 1 near the arc {( = e ia ,a G [0, |vr]}, 

(ib) 0o(C) ^ near the complementary arc {( = e ia ,a G [^tt, 2tt}}; 

(ii) 0o is subharmonic in B(0, |), moreover [A0 o ](C) ^ const > for all ( G -8(0, |); 

(iii) O (C) ^ const > for all C G B(0, §). 

Clearly, O exists, if the constants in (ii),(iii) are small enough. Let 

tf( u ) : =J^A for ueT . 

\a - w y 

Then, (f> 5 ^ 1 on the discrete arc afe and (p s ^ on the complementary arc 6a. 

If 5/r is small enough, then, due to (ii) and the approximation property (Lemma l2.2l) . 
cj) 5 is discrete subharmonic in -Bf(wo,r). Using the maximum principle, one obtains 

u s (u; ab; B s v (0, r)) ^ <p s {u) ^ const > for all u G B*(0, \r). 

If S/r ^ const > 0, then the claim is trivial, since the random walk starting at uq can 
reach the discrete arc ab in a uniformly bounded number of steps. □ 



cxp 



— const ■ - 



R-r 



^ 777 — \ < ex P 
H(ui) 



10 



DMITRY CHELKAK AND STANISLAV SMIRNOV 




Figure 3. (A) A discrete disc. The "black" polygonal boundary B 
and the "white" contour W are shown together with the correspondences 
z i — > b(z), z G W PI <0, and z \— > w(z), z G B fl <^>. (B) The proof of the 
weak Beurling-type estimate (Proposition 12.101) . The probability that 
the random walk makes a whole turn inside the annulus (and so hits the 
boundary dfl 6 ) is uniformly bounded away from due to Lemma [2.91 



Let f2p be some connected discrete domain, u G f2f and E C dVt 6 T . We set 



dist^ (u; E) = inf{i? : u and E are connected in f2 r f] B(u } R)}. 

The following Proposition is a discrete version of the classical Beurling estimate with 
power 1/2 replaced by some positive (3. 

Proposition 2.10 (weak Beurling-type estimates). There exists an absolute con- 
stant (3 > such that for any simply connected discrete domain flf^, point u G Int fif 
and some part of the boundary E C <9fi r one has 



u s (u;E;tt s r ) ^ const 



dist(w;<9fif) 



P 



_dist n «(u;£) 

Above we set di&mE := 5, if E consists of one point 



and uj 5 (u;E;Q 5 r ) ^ const- 



di&mE 
_dist n s (u;E) 



Proof. The proof is quite standard. We prove the first estimate, the second is similar. 
Let d = dist(w; <9f2 r ) and r = dist c « (u; E). Without loss of generality, one may assume 

that C ■ 5 ^ d ^ |r for a sufficiently large constant C . The quantity E] fif) 

is the probability that the random walk starting at u first hits the boundary of Jl r 
inside E. Using Lemma [2.91 (see Fig. [33), it is easy to show that for each d ^ r' ^ |r 
the probability to cross the annulus B{u, 2r') \ B(u, r') inside f2 r without touching its 
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boundary is bounded above by some absolute constant p < 1 that does not depend on 
r' and the shape of f2f . Hence, 

u\u-E-n 5 r ) <: p l °S2(r/d)-l =p -i . (rf/ r )-i°g 2 p ; 
so one can take the exponent /3 = — log 2 p > 0. □ 

2.4. Factorization of the discrete Laplacian. Discrete holomorphic functions. 

Let <) denotes the set of the rhombi centers. Above we discussed the theory of discrete 
harmonic functions defined on the isoradial graph T (or, in a similar manner, on its 
dual r*). Now, following [MerOlj and |Ken02j . we introduce the notion of discrete 
holomorphic functions. These are defined either on A = T U T*, or on <^ = A*. Note 
that, in contrast to similar T and T*, A and <0 have essentially different combinatorial 
properties. 

Definition 2.11. Let z E (} be the center of the rhombus V1V2V3V4, where the vertices 
fi,t>3 G T, f2,t>4 £ T* ore enumerated in counter clockwise order. We define the 
(antisymmetric) weights of the edges {vjz), (zVj) by 

Hv,z = i(v j+1 -Vj-i), /j, zv . = i(vj-i-Vj + x) 

and the weight of z G (} (the center of the rhombus with half-angles 0(z), |7T — 9{z)) by 

//I (2) = Area(t>if 2 f3f4) = S 2 sm26(z). 

Foru G T (and, in the same way, forw G T*) we set fi A (u) = |yUp(u) = | XL S ~ U /-4>( z s)- 

Clearly, formulas similar to (12.21) are fulfilled for defined on subsets of ^> or A. 

Definition 2.12. (i) Let F be defined on (some part of) (}. We introduce its discrete 
derivatives d 5 F and d 5 F (defined at v G A) as 

[d s F](v) = -j— \2^ z F(z s ) and [d s F](v) = — * £ J^F(z s ). 



We call F discrete holomorphic at v iff [d s F](v) = 0. 

(ii) Let H be defined on (some part of) A. We define d s H and d s H (at z G (}) as 



v j r ^ z 



H( Vl )-H(vz) | H(v 2 )-H(v 4 ) 

V1—V3 V2-V4 
H( Vl )-H(v 3 ) | H(v 2 )-H(v,) 



We call H discrete holomorphic at z iff [d s H](z) = 0. We use the same notations, 
if H is defined on T (or T*) only. In this case we formally set H\r* = (or H\r = 0). 

This definition is the natural discretization of the identities 
(dF)(x+iy)dxdy = -- I F(()d( and / / (dF)(x+iy)dxdy = % - <f F(()d(, 

1 JdP J J P 1 JdP 
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where P denotes the dual face W\W 2 ..w n containing v or the rectangle v \V 2 v 31*4, where 
Vj = Vj + Vj + \ — z (note that the area of P is /if( u ) — ^^a( v ) or 2fJ>%(z), respectively). 
Thus, d s and d 5 have approximation properties similar to those in Lemma l2~2l Namely, 

W<l>\<>](v)- |[S*0W(«)-(S0)(t;)| = 0(5), u e A, 

PV|a](^)-(90)(^)|, |PVIa](*)-(8M*)| = 0(5 2 ), z£ 

for smooth functions (the additional cancellation for d s (4>\\) and c^(</>|a) comes from 
the symmetry of the rectangle and, in general, does not hold for d s ((j)\^), d s ((p\^)). The 
following factorization of the discrete Laplacian was noted in jMerOl] and [Ken02| . 

Proposition 2.13. For functions H defined on subsets of A the following is fulfilled: 

[A s H](u) = A[d s d s H]{u) = A[d 5 d s H](u) 

at all points u £ A where the right-hand sides make sense. 

Proof. Straightforward computations show (see Fig. [IB) 

[d s d s H}(u) = — * [t^6 s ■ [H(u s )-H(u)] - i ■ [H(w s+1 )-H(w s )]] = I^SM 

and similarly for [d s d s H](u). □ 

2.5. Elementary properties of discrete holomorphic functions. Let H and F 

be defined on some subsets of Y (or T* or A) and ^>, respectively. 

(1) If H is holomorphic on A, then H is harmonic, i.e. H\r and H\r* are harmonic. 

(2) Conversely, in simply connected domains, H is harmonic on T if and only if 
there exists a harmonic function H on T* such that H + iH (as a function on 
A) is holomorphic. H is defined uniquely up to an additive constant. 

(3) If H is harmonic on A, then its derivative F = d s H is holomorphic on 

(4) In simply connected domains, if F is holomorphic on <0>, then there exists a 
function H = J 5 F(z)d s z harmonic on A such that d 5 H = F. Its components 
H\-p and H\r* are defined uniquely up to additive constants by 

H(u 2 ) - H(ui) = F (~(u 2 +ui)) ■ (u 2 - u 2 ~ Mi, 

where U\, u 2 £ Y or u\, u 2 £ T*, respectively. 

(5) F is holomorphic on ^> if and only if both functions 



[BF)(z) := Pr ^F(z);ui(z)—u 2 (z) J , [WF)(z) :—Px F(z); w 1 (z) — w 2 (z) 

are holomorphic, where Ui j2 (z) £ T and Wi t2 (z) £ T* are the black and white 
neighbors of z £ <C>, respectively, and Pr[-;^] denotes the orthogonal projection 
onto the line £1R (note that F = BF + WF, since u\(z)— u 2 {z) J_ W\(z)—w 2 {z)) . 

(6) If H is real and harmonic on T (or T*) and F = d s H (recall that F is holomor- 
phic on <^), then F = BF (or F = WF, respectively). 

(7) If H is holomorphic on A, then the averaged function 

[m 5 H)(z) := - [H( Ul )+H(u2)+H(u 3 )+H(ut)], z ~ Uj £ A (2.6) 
is holomorphic on {>. 
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Proof. (1) Easily follows by writing A S H = 4d 5 d s H = 0. 

(2) [A s H](u) = iff the sum of increments H(w s+ i) — H(w s ) around u is zero. 

(3) Easily follows by writing d s F = ^B s d 5 H = 0. 

(4) H is well-defined due to d 5 F = 0. By definition, d s H = F, so H is harmonic. 

(5) Clear, since Re\d 6 F] = d 5 [BF] and lm[d s F] = d 5 [WF}. 

(6) Directly follows from definitions. 

(7) Straightforward calculations show [d s m 5 H]{u) = 0, since [d 5 H](z s ) = implies 



Below we will also need the averaging operator m 6 (adjoint to (12. 6p ) for functions 
defined on (some part of) <^: 



Unfortunately, there are several unpleasant facts that make discrete complex analysis 
on rhombic lattices more complicated than the standard continuous theory and even 
than the square lattice discretization: 

• One cannot (pointwise) multiply discrete holomorphic functions: the product 
FG is not necessary holomorphic if both F and G are holomorphic. 

• One cannot differentiate discrete holomorphic functions infinitely many times: 
d^F is not necessary holomorphic (on A) if F is holomorphic on <^. Similarly, 
d 5 d 5 H is not necessary holomorphic (on A) if H is holomorphic on A. 

• There are no "local" discrete analogues of d that map holomorphic functions 
on A (or on ^>) to holomorphic functions defined on the same set (A or ^>). One 
cannot use the natural combination of m 5 and d 5 since, in general, d 5 does not 
map Hol(<)) into Hol(A) and m s does not map Hol(A) into Hol(^>). 

The first obstacle (multiplication) exists in all discrete theories. Concerning the sec- 
ond, note that in our case there is some "nonlocal" discrete differentiation (so-called 
dual integration, see |Duf6 8] and |Mer07j ). Also in two particular cases the local differ- 
entiation leads to holomorphic function again: for the classical definition on the square 
grid (see the book by J. Lelong-Ferrand [L-F55j ) and for a particular definition on the 
triangular lattice [DN03]. 

2.6. The Cauchy kernel. The Cauchy formula. Lipschitzness. The following 
asymptotic form of the discrete Cauchy kernel is due to R. Kenyon. 

Theorem 2.14 (Kenyon). Let z G <^. There exists a unique function F = K(-; z ) : 
A — >• C such that 

(i) [d s F}(z) = for all z ^ z and ^(z ) • {d s F](z ) = 1; 

(ii) \F(u)\ — > as \u— Zq\ — > oo. 
Moreover, the following asymptotics holds: 



[m s H}(z s ) 



H{u) H{w s+ i)(w s+ i-u) - H{w s ){w s -u) 
2 2{w s+1 -w s ) 



□ 




(2.7) 
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K(w] Zq) = — Pr 



7T 



(!) ...(2) 



) w 







+ 



5 



w g r, 



W-Zq 

where u^\u^ G T and Wq,Wq G T* are the black and white neighbors of z , respec- 
tively. 

Proof. We give a short sketch of R. Kenyon's arguments [Ken02j (together with a slight 
but necessary modification) in Appendix IA.1I □ 

Let fif be a bounded simply connected discrete domain (see Fig. [2A,[3]A.). Denote 
by B = UQUiU2-.u n , u s G T, its closed polyline boundary, enumerated in counter 
clockwise order. Denote by W = WoWxW2--W m , w s G T*, the closed polyline path 
passing throw the centers of all faces touching B (enumerated in counter clockwise 
order). For functions G defined onBn<) (or, in the same way, on W H (}) we set 



<t G{z)d 5 z := VG(|( Ms+ i+tt s )) • K+i- 



be a discrete 



Lemma 2.15 (Cauchy formula). Let tt 6 ^ = tt 5 n <>, F : fi^ -> 

holomorphic function (i.e., (d F)(v) = for all v e Q 5 H A) and zq G fi^ \(BU W). 
Tnen 



*■(*> = Ti 



B 



K(w{z)]z )F(z)d 8 z+ <f> K(u{z)]z Q )F(z)d s z 

Jw 

where z ~ w(z) G W , if z G B n 0, and 2 ~ G5, i/z G W n fsee Fza. 
Proof. Discrete integration by parts and [d s F](z) = give 



neBUVK 



7 K(v;z )[i vz F(z) +j ^ K(v;z )fi vz F(z). 



v=u(z), z£W 



v=w(z), z£B 



It's easy to see that n u t z ) z = i(w s+ i — w s ) for all z = |(w s +i + w s ) G W fl <0> and 



fJ>w(z)z 



-l[Us 



-i— u s ) for all z = |(m s+ i+m s ) G B n <>. 



□ 



If F — BF, i.e., F is the discrete gradient of some discrete harmonic function 
H : fip — * R (or, in a similar manner, if F — WF), then the Cauchy formula may be 
nicely rewritten in the asymptotic form: 

Corollary 2.16. Let F = BF : -> C 5e discrete holomorphic and z G ft£\(.BUW). 
Tnen 



F(^o) = Pr 



1 

2m 



F(z)d 5 i 



+ 



F(z)d 5 ; 



; «o 



(!) ..(2) 



a 



o 



+ 



(M5L\ 

WW 



' B 2 ~ ^0 JVF 2 — ^0 

where d = dist(zo, W), M = max ze Buw \F(z)\ and L is the length of B U W. 
Proof. We plug Kenyon's asymptotics (Theorem 12 ,14ft of K(v;zq) into Lemma [2. 151 
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If z G B n 0, then F(z)d s z G R, so 



Zq) — = Pr 



F (z)d S z . (1) (2) 



, u u 



O 



since Pr[/; r] • £ = Pr[/£; ir] for £ G zM and Wq 1 "*— Wq 2 ' 1 -L UQ — Wo 



27fi(z — Zq) 

' "u — "u 



If z G IV n 0, then F(z) ■ cf z G z'R and 



zq) — — = Pr 



F(z)d s z _ A1) A2) 



2m(z-z ) ' U ° U ° 



+ 



M5\d s : 
dP~ 



M5\d s ; 



d 2 

since Pr[/; t] • £ = Pr[/£;r] for £ G R. Summarizing, one obtains the result. □ 

Proposition 2.17 (Lipschitzness of discrete holomorphic functions). Letuo G T 
and let F = BF be discrete holomorphic in B^(uq, R). Then 

\F(z s ) — Pr \2[m s F](uo); u s — u ] I ^ const ■— -, where M— max l-Ff^ll, 

" ' R B^u^R) 

for all z s = ~(u s +Uo) ~ Uq. Furthermore, ifu±,u 2 G B^(uQ,r) C Int By(uo, R), then 
\[m 5 F]{u 2 ) - [m 5 F}( Ul )\ < const M ' M2-Ml ' 



R-r 

where the averaging operator m 5 is defined by p£ 



Proof. Let B and W be the same discrete contours as above (see Fig. [3K). Note 
that their lengths are bounded by const -R. Applying the asymptotic Cauchy formula 
(Corollary 12.161) for all z s ~ u , one obtains 

„, \ ^ rA . ^{M5\ , , 1 f f F(z)d s z f F(z)d s z 

F(z s ) = Pr[A;u s -u ] +0[ — ), where A = —[<p K J - + 1 

Due to the identity 



R J 27U \Jb Z — Uq J w Z — Uq 



^(2.)Pr[A;«.-« ] = — jy-r 5 2 sin2^ 



4<W^ v L J 4 <H 

2 + 16z/il(«) ^ [ ' 2 



(see Fig. [TB), it gives 

[-^] W ^ + (f 

In particular, |-F(z a ) — Pi[2[m s F](u ); u s — Uq\ | ^ const -M5/R. Using similar calcula- 
tions for \m 5 F]{u s ) and the trivial estimate \(z— u s ) _1 — (2— m ) _1 | ^ const -S/R 2 , one 
obtains 

|[m 5 F](M s ) — [m <5 F](M )| ^ const -—5-, if ~ ^o- 

Clearly, the estimate for | [m s F](u 2 ) — [m 5 F](ui) \ follows by summation along the path 
(with length ^ const -\u2— v>i\/5, which exists by condition (4)) connecting U\ and u 2 
inside B^(u , r). □ 
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3. Convergence theorems 

3.1. Precompactness in the C 1 — topology. In the continuous setup, each uniformly 
bounded family of harmonic functions (defined in some common domain fl) is precom- 
pact in the C°°-topology. Using Corollary 12.81 and Proposition 12 .171 it is easy to prove 
the analogue of this statement for discrete harmonic functions. Let H s : f2p — > R be 
defined in f2p C T s . Then, as usual, H s can be thought of as defined in Q 5 C C by 
some standard continuation (say, linear on edges and harmonic inside faces). 

Proposition 3.1. Let H Sj : f2p — > R be discrete harmonic functions defined in discrete 
domains with 5j^0. Let C (JS f^Z ^ c C be 

some continuous domain. If 

are uniformly bounded on Q, i.e. 

max *, \H Sj (u)\ < M < +oo for all j, 

then there exists a subsequence 5j k — > (which we denote by 5k for shortness) and two 
functions h : O — > R ; / : Q — > C suc/i i/ia£ (u>e denote by uniform convergence) 

H 5k =4 uniformly on compact subsets K C 



and 



i t Fi 

|«2 — «i| 



Re 



/(«) 



{• k 

U 2 — «1 

1 k k 

\U 2 — ^1 | J 



(3.1) 



if u\, u\ G T 5fc , u\ ~ and w^, u\ — > w G K C fi. Moreover, the limit function h is 
harmonic in Q, \h\ ^ M and f — h' x — ih' y = 2dh is analytic in Q. 

Remark 3.2. In other words, the discrete gradients of H s defined by the left-hand side 
of Ii3. 1\) converge to V/i. Looking at the edge {u\u 2 ) one sees only the discrete directional 
derivative of H 5 along r = (u 2 — Ui)/\u 2 — ui\ which converges to (Vh(u),r) = 2dh{u)-r. 

Proof. Due to the uniform Lipschitzness of bounded discrete harmonic functions (see 
Corollary ES]) and the Arzela-Ascoli Theorem, the sequence {H Sj } is precompact in the 
uniform topology on any compact subset K C Vt. Moreover, their discrete derivatives 
(defined for z G Q S £) 

are discrete holomorphic and bounded on any compact subset K C Q. Then, due to 
Proposition l2.17l (ii) and the Arzela-Ascoli Theorem, the sequence of averaged functions 
m Sj F Sj (defined on fip by (12.71) ) is precompact in the uniform topology on any compact 
subset of Q. Thus, for some subsequence 8k — > 0, one has 

H Sk =4 h and 2m 5k F &k =4 / 

uniformly on compact subsets of Q. Moreover, due to Proposition 12 . 1 71 (i) . it gives 

F Sk (z) - Pr 



f(z);u 2 (z) -ui(z) 



uniformly on compact subsets of f2. 



It's easy to see that h is harmonic. Indeed, let <fi : Q — > R be an arbitrary C^°(Q) test 
function (i.e., G C°° and supp </> C fi). Denote by h s , (p s and (A0) 5 the restrictions of 
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h, 4> and Acf> onto the lattice V s . The approximation properties ( (12.21) and Lemma 12.21) 
and discrete integration by parts give 

(h,A4>) a = lim V^H(A0) 5 H^(n) = lk V^WWW 

= Jg% E # 5 («)[AV](«)/4(«) = Jim Q ^ [A 5 tf a ]( U )0V)/4(«) = 0. 
Furthermore, for any path [MoWn] 5 = u oUi-.u n , u s+ ± ~ u s , u s G T 5 , one has 

/■J n-l 

#(« n ) - ff(«o) = / F%z)d 5 z = ^F j (i(« s+1 +« s )) ■ (n s+1 -n s ). 

Taking appropriate discrete approximations of segments [w] C Q (recall that rhombi 
angles are bounded from and tt, so one may find polyline approximations with uni- 
formly bounded lengths) and passing to the limit as 5 = 5k — > 0, one obtains 



h(v) - h(u) = I Re[f{z)dz) = Re / 

J \uv\ J \% 



f{z)dz 



for all segments [uv] C f2. 



It gives ah' x (u)+(3h' y (u) = Re[(a+if3)f(u)] for all u G Q and a, /3 G K, so / = 29/i. □ 

3.2. Caratheodory topology and uniform C 1 -convergence. Below we need some 
standard concepts of geometric function theory (see [Pom 92] . Chapters 1,2). 

Let Q be a simply connected domain. A crosscut C of Q is an open Jordan arc in 
Q such that C = C U {a, b} with a, b G <9f2. A prime end of f2 is an equivalence class of 
sequences (null-chains) (C n ) such that C n nC„ + i = 0, C n separates Co from C7 n+ i and 
diamC n — > as n — > oo (two null chains (C n ), (C)„ are equivalent iff for all sufficiently 
large m there exists n such that C m separates Co from C n and C m separates Co from 

C n ). 

Let -P(^) denote the set of all prime ends of Q and let <fi : Q — > D be a conformal 
map. Then (see Theorem 2.15 in [Pom92j) <p induces the natural bijection between 
P(O) and the unit circle T = <9D. 

Let Mo £ C be given and Q n , Q C C, be simply connected domains ^ C with 
Uo £ ^- We say that Sl„ ^ O as n ->oo m i/ie sense o/ kernel convergence with 
respect to u iff 

(i) some neighborhood of every u G f2 lies in fi n for large enough n; 

(ii) for every a G <9f2 there exist a n G <9fi n such that a n — > a as n — > oo. + 

Let 0fc : fifc — ► D, : — > D be the Riemann uniformization maps normalized at u 
(i.e., <p(uo) = and 4>'(u ) > 0). Then (see Theorem 1.8 [Pom92j) 

ilfc — » w.r.t. «o -v^ ^^T 1 =4 _1 uniformly on compact subsets of D. 

It is easy to see that (a) <pk =4 as — > oo uniformly on compact subsets if C Q; 
(b) for all < r < i? the set of simply connected domains Q : B{uq, r) C Q C 5(0, i?) 
is compact in the topology of kernel convergence w.r.t. uq. 

Definition 3.3. Let Q = (fi; v, ..; a, b, ..) be a simply connected bounded domain with 
several (possibly none) marked interior points v, .. G Int Q and prime ends (boundary 
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points) a,b, .. G -P(fl) (we admit coincident marked points, say, a = b) and let u G fl. 
We write 

(fl fc ; u k ) = (fit; u k , v k , ••; a k , b k , ..) (Q;u) = (Q;u,v, ..;a,b, ..) as fc — > oo, 

ijf £/ie domains Q k are uniformly bounded, — > w ; fl^ — > fl in the sense of kernel 
convergence w.r.t. u and 4>h{ v k) —> 4>{ v )i--> 0fc( a fc) ~^ 4>{ a )i--> where 4>k : fl& — > D, 
: fl — > D are i/ie Riemann uniformization maps normalized at u. 

Remark 3.4. Since one /ms 1 0(f) I < 1. T/ras, (f>k{v k ) — > 0(f) implies Vk — > v. 

Moreover, one can equivalently use the point v instead of u in the definition given above. 

Definition 3.5. Let fl be a simply connected bounded domain, u,v, .. G fl and r > 0. 
We say that the inner points u, v, .. are jointly r-inside fl iff B(u, r), B(v, r), .. C fl 
and there are paths L uv , .. connecting these points r-inside fl (i.e., dist(L uv , <9fl), .. ^ r). 
In other words, u,v, .. belong to the same connected component of the r-interior o/fl. 

Note that for each < r < R there exists some C(r, R) such that, if fl C -8(0, R) 
and u, v, .. are jointly r-inside fl, then 

\<j>(v)\,..^C(r,R)<l, (3.2) 

where : fl — > © is the Riemann uniformization map normalized at u. Indeed, con- 
sidering the standard plane metric, one concludes that the extremal distance between 
L uv and <9fl in fl \ L uv is no greater than irR 2 /r. Thus, the conformal modulus of the 
annulus D \ <p(L uv ) is bounded below by some const (r, R) > 0. Since <p(u) = 0, (13.21) 
holds true for some C(r,R) < 1. 

Suppose that a harmonic function 

h( - ; fl) = h( • , v , ..; a, 6, ..; fl) : fl -»• R 

is defined for each fl = (fl; v, ..; a, b, ..). Similarly, let flf = (fl*; f 5 , ..; a 5 , o 5 , ..) denote a 
simply connected bounded discrete domain with several marked vertices v s , .. G Int flf 
and a s ,b 5 ,.. G <9flf and 

H 5 ( • ; flp) = # 5 ( • , v s , ..; a s , b 5 , ..; flf) : flf -> R 

be some discrete harmonic function in flp. Recall that fl 5 denotes the corresponding 
continuous polygonal domain. 

Definition 3.6. We say that the functions H 5 are uniformly C 1 — close to h inside 
fl 5 , iff for all < r < R there exists e(5) = e(5,r,R) — > as 5 — > such that for all 
discrete domains flf = (fif; v s , ..; a s , b s , ..) and u 5 G Int flp the following holds true: 

//fl 5 C B(0,R) andu s ,v s ,.. are jointly r-inside VI s , then 

\H s (u s ,v s ,..;a s } b 5 ,..;n 5 r ) - h(u s , v 5 , ..; a 5 , b s , ..; fl 5 ) | ^ e(8) (3.3) 
and, for all u 8 ~ u\ G fif, 



iJ 5 (uf ; flf) - # 5 (f/; flf) 



If 5 - u 8 \ 



where 2dh = h' x — ih' y . 



Re 



2dh(u s ;Q s 



u{ - u 8 
\v$ - it 5 1 



< e(S), (3.4) 
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Proposition 3.7. Let (a) h be Caratheodory-stable, i.e., 

h(uk',Qk) ~~ * h(u;Q), if (ft^w/J — 5 (ft;w) as k — > oo; (3-5) 
and (b) H 5 — > h pointwise as 5 — > 0, i.e., 

H s (u s ;n 5 r ) ->h(u;Sl), if (ft 5 ; u 5 ) ^5 (ft; m) as 5^0. (3.6) 

TTien i/ie functions H 5 are uniformly C 1 -close to h inside ft 5 (see Definition ^. 6\) . 

Proof. Assume ( 13.31) does not hold true, i.e., 

\H S (u s ; ftf) - /i(w 5 ;ft 5 )| ^e >0 

for some sequence (ftf; u s ), 5 = 5j —> 0, such that B(u s , r) C ft 5 C -8(0, R). Taking a 
subsequence, one may assume that u s — > u for some w G B(0,R). The set of all simply 
connected domains ft : B(u, |r) C ft C -8(0, R) is compact in the Caratheodory 
topology. Thus, taking a subsequence again, one may assume that 

(ft 5 ; u 6 , v 5 , ..; a 5 , b s , ..) ^5 (ft; M , u, ..; a, 6, ..) as 5 = 5 fe -> 

(note that the marked points v 5 ,.. cannot reach the boundary due to (13.21) ). Then, 
(a) the Caratheodory-stability of h and (b) the pointwise convergence H s — > h give a 
contradiction. Indeed, both 

(a) h(u 5 ; ft 5 ) -> %; ft) and (b) F 5 (u 5 ; ftf) -»• %; ft) as 6 = 5 k -> 0. 

In view of Proposition 13.11 the proof for discrete gradients goes by the same way. 
Assume ( 13 ,4p does not hold for some sequence of discrete domains. As above, one may 

take a subsequence 5 = 5 k such that (ft 5 ; u 5 ) -^5(ft, u). The Caratheodory-stability 
implies 

h( ■ ; ft 5 ) =4 h( ■ ; ft) uniformly on ~r) as 5 = — > 0. 

Indeed, (ft 5 ; 5) ^5(ft; 5) for all u G 5(w, |r). If |/i(w 5 ; ft 5 ) -/i(w 5 ; ft) | ^ e > for some 
u s and all 5 = 5 k , then, taking a subsequence 5 = 5 m such that u s — > 2 G -B(w, |r), one 
obtains a contradiction with /i(-u 5 ; ft 5 ) — > ^(u; ft) and /i(m 5 ; ft) — > h(u; ft). 

The uniform estimate \H 5 ( ■ ; ftf) - /i( • ; ft 5 )| ^ e(<J) -> (see above) gives 

# 5 ( • ; ftf) =4 /i( • ; ft) uniformly on B{u,\r) as 5 = 5 k -> 0. 

In particular, the functions if 5fe (-; ftf fc ) are uniformly bounded in -B(it, |r). Thus, using 
Proposition 13.11 one can find a subsequence 5 = 5 m such that the discrete derivatives 
of H Sm converge (as defined by ( 13.41) ) to / = 2dh(-; ft) which gives a contradiction. □ 

3.3. Convergence for the harmonic measure of arcs and the Green's function. 

Theorem 3.8. The discrete harmonic measures o/(-; b s a 5 ; ftp) are uniformly C x -close 
to the continuous harmonic measures uj{-; a s b s ; ft 5 ) inside ft 5 (see Definition ^. 6}) . 

Proof. By conformal invariance, the continuous harmonic measure is Caratheodory 
stable, so the first assumption in Proposition 13.71 holds true. Thus, it is sufficient to 

prove pointwise convergence ( \3.6h . Let (ft 5 ;w 5 ;a 5 , b 5 ) ^5(ft ;u;a,b). Obviously, the 
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functions u/(-; b s a 5 ; f2p) are uniformly bounded in Q. Due to Proposition 13.11 one can 
find a subsequence 5k — > such that 

u Sk (-,b Sk a Sk ;n 5 r h ) =4 H 

uniformly on compact subsets of Q, where H : Q — > R is some harmonic function. It 
is sufficient to prove that H{u) = u(u; ba; fi) for each subsequential limit. 

Let z s G IntOp. The weak Beurling-type estimate (see Proposition 12.101 ) gives 

dist(z s ;dtt s r ) 1/3 



< u 5 (z 5 ; b 5 a 5 ; Q 5 r ) < const ■ 



dist n « (z s ; b 5 a s 



uniformly as 5 — > 0. Passing to the limit as 5 = 5k — > 0, one obtains 

dist(z;dtt) lP 



< H(z) < const 



for all z G Int f2. 



distn(z; 6a) 

Therefore, H = on the arc a6 C <9fi. Similar arguments give H = 1 on the arc 
6a C dfl. Hence, H = u>(-; ba; fl) and, in particular, H(u) = u{u; ba; ft). □ 

Let fif be some bounded simply connected discrete domain. Recall that the 
discrete Green's function ■ ; v s ), v s G Int fi r , is given by 

G n s r ( • y ) = G r ( • ; - G* * ( ■ ;v 5 ) : fi r -> R 

(see (12.41) ). Theorem 12.51 claims uniform C 1 -convergence of the free Green's function 
Gr to its continuous counterpart 

G c (u;v) = -!- log |u-v| 

with an error of order 0(5 2 \u — v\~ 2 ) for the functions and so 0(5\u — v |~ 2 ) for their 
gradients. Let Gq(-; v) = Gc{-', v) — Gq(-; v) : Q — > R be the solution of the Dirichlet 
problem 

AG* n (-;v) = in ft and G* Q (-; v) = G c (-, v) on ffl. 

Theorem 3.9. The discrete harmonic functions G^ s (-',v s ) are uniformly C 1 -close to 
their continuous counterparts G^-; v 5 ) inside Q s (see Definition ^. 6\) . 

Proof. We start with the proof of lETE|l . Let as 5 = ^ -> 0. 

We need to check that G^,, (m 5 ; i/) — > G^(u;^). Recall that — > u, v s — > t> and 

-B(w, r), -B(f , r) C ft 5 for some r > 0, if 5 is small enough. Let 

g(z) = Gc(z;v) for z : |z— v| ^ \r. 

Denote by £?*(•; Q r ) the solution of the discrete Dirichlet problem 

A 5 G* = in n s r and G* = g on <9ft£. 

Due to the maximum principle, since v 5 — > v and \Gr(-;v 5 ) — Gc(-] v 5 )\ — ► as 5 — ► 
uniformly outside -B(v, |r), it is sufficient to prove that 

G*(w 5 ; flf ) -> Q(w; u) as 5 = 5j -> 0. 
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Note that the functions G*(-;fif) are uniformly bounded. It follows from Proposi- 
tion 13.11 that one can find some subsequence 5k — > such that 

G*{-;Q 5 r k ) =4 H 

uniformly on compact subsets of Q, where H : Q — > R is some harmonic function. It 
is sufficient to prove that H{u) = Gq(u; v) for each subsequential limit H. 

Let z s G fif flint Q, \z s —v\ ^ r and dist(z 5 ; dfl) = d. For each d < p < R one has 

G*{z s ;Q s r )-g(z s ) < max -#(><% w a e n B(z s , p)} 

+ max{|^(w' 5 )-g(z <5 )|, G <9ftf} • w*^; dfif \ £(/, p); ftf ). 

Since f2f are uniformly bounded by our assumption (see Definition 13.61) . the function 
g is uniformly bounded and Lipschitz in Q s \ B(v, |r). Moreover, dist^ 5 ; <9f2f) ^ 2d, 
if 5 is small enough. Applying Proposition 12.101 one obtains 

G*(z s ; tt s r )-g(z s ) = 0(p + {2d/pf) uniformly as 8 = 5 k -»• 0. 

Setting p = df ) + 1 and passing to the limit as 5 = 5 k — > 0, one obtains 

jtf^) -#(z)| = O ((dist(z,dft))/&) for all z G \ \r). 

Thus, H = g = Gq(-; v) on the boundary dQ. Since both H and Gq(-; v) are harmonic 
in Q, it gives i?(tt) = Gq(u; v). 

Repeating the arguments given above in the continuous setup, it's easy to derive 
the Caratheodory stability of the function Gq(-;v). Thus, both assumptions (13.51) and 
in Proposition 13.71 hold true. □ 



3.4. Convergence for the Poisson kernel normalized at an interior point. Let 

f2 r be some simply connected discrete domain, a s G <9f2 r and v s G Int f2 r . We call 

P s = P 5 (-;v 5 ;a s ;n s r ) : tt 5 r 
the discrete Poisson kernel normalized at v s if 



<5o<5 _ n :„ o<5 o<5| q and nSf.M 

LU s (v s ;{a s };n 8 r ) 



A'P* = 0innJ, n fl o f \ { a} = ° ™ (1 p? f 

Note that 



V s ; a s ; fif) 



In the continuous setup, let Q be some simply connected domain, a G P(fi) be 
some prime end and v G Int Q. We call P = P(-; v; a; O) the Poisson kernel normalized 
at v, if 



AP = in fi, P lan\{a}= ' P >° and ^0) = 1 

(note that P is uniquely defined by these conditions for arbitrary simply connected 
domains Q as the conformal image of the standard Poisson kernel defined in D). 

Theorem 3.10. The discrete Poisson kernels P s (-; v 5 ; a 5 ; fi r ) are uniformly C 1 -close 
to the continuous Poisson kernels P(-; v ; a 5 ; VI s ) inside VI s (see Definition ^. 6]) . 
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Figure 4. Parts of the continuous domain fl and some discrete domain 
fl s close to fl in the Caratheodory topology. Marked boundary points 
a s G dfl s and a G P(fl) are close to each other in this topology since 
there are corresponding small cross-cuts near a d . L d denotes the closest 
to v arc in {z : \z — a d \ = d} fl fl 5 . The quadrilateral R d d is shaded. 

Proof. The continuous Poisson kernel is Caratheodory stable by definition, so ( 13.51) 
holds true. Thus, it is sufficient to prove pointwise convergence (13.61) . 

Let (fl 5 ; u 5 , v s ; a 5 ) ^-5(£7;w, v;a). Recall that v 5 — > v and B(v,r) C fl 5 for some 
r > 0, if 5 is small enough. It follows from P(v s ) = 1 and the discrete Harnack Lemma 
(Proposition 12.71 (ii)) that P s are uniformly bounded on each compact subset of fl. 
Then, due to Proposition 13.11 one can find a subsequence 5k — > such that 

P 5k (-;v 5k ;a 5k ;tt 5 r k ) =4 H 

uniformly on compact subsets of Q, where H ^ is some harmonic function in Q. It 
is sufficient to prove that H(u) = P(u\ v; a; fl) for each subsequential limit H. 

Let d > be small enough. Then, there exists a crosscut 7^ C B(a d , |d) in fl 
separating a and u,v (see Fig. HI). Moreover, one may assume that u, v £ B{a d ^d) 
and u, v belong to the same component of fl \ B(a d , Ad). For sufficiently small 5 let 

L d C {z : \z-a d \ = d} H fl 5 

be an arc separating v s and a s in fl 5 (we take the arc closest to v 5 , see Fig. HI). Let fl 5 d 

denote the connected component of fl 5 \ L d containing v 5 . Since fl s ^5 fl w.r.t. v and 
v s — > v , one has 

u(v 5 ; L d ; fl 5 d ) ^ const(rf) > (3.7) 

(here and below constants const (d) do not depend on S). Similarly, let fl^ d be the 
connected component of fl s \ L^ d containing u, v. Denote 

M 5 3d = m a x{P 5 (z 5 ), z 5 efi 6 M nv 6 }. 
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Since the function P 5 is discrete harmonic, one has 

for some nearest- neig hbor path K S M = {z$ ~ z{ ~ z 5 2 ~ .., z s s G fif }, starting at some 
z o e ^ld- Since -P <5 |an*\{a< 5 } = 0; the unique possibility for this path to end is a 5 . 

Using (13. 7p . it is not hard to conclude (see Lemma [5.111 below) that the following 
holds true for the continuous harmonic measures: 

lu(v s ; K 3d ; tt 5 \ K 3d ) ^ lu(v s ; K 3d H Q s d \ K 3d ) ^ const -lu(v 5 ; L d ; n s d ) > const(d), 

where K 3d is the corresponding polyline starting at Zq and ending at a s . In view of 
Theorem 13.81 if 5 is small enough, then the same inequality 

u s (v s ; K s 3d ; Q s r \ K* u ) > const(rf) > 

holds true for discrete harmonic measures uniformly as 5 — > (with smaller const {d)). 
Recall that P 5 (v s ) = 1 by definition and P s (v) ^ M 5 U along the path Mf d . Thus, 

Mg d ^ const (ci), if 5 is small enough. 

Finally, let z s G n f2g d be such that l^" 5 — a d \ > 3d. The weak Beurling-type estimate 
immediately gives 



P 5 (z 5 ) ^ const 



dist(z 5 ;dQ 5 ) 



_dist{z s ;L 3d ) _ 
Passing to the limit as 5 = 5k — > 0, one obtains 



Ml d ^ const (d) 



d\st{z 5 ]dQ 5 )' 
1 — Clrf I — 3d 



H{z) s= const (d) 



dist(z; <9Q) 



T/3 



for all z G £l 3d suchthat |z— aJ > 3d. 



\z — a d \ — 3d 

Thus, H = on dVl \ {a}. Since H > and /f(v) = l, it gives H = P(-;v, a; fi). □ 

Lemma 3.11. Lei C C 6e some simply connected domain, v G Int Q and a G P(fi). 
Let L d C {z : \z—a d \ = d}C\Q be the arc separating v and a that is closest to v, and Q d 
be the connected component offl\L d containing v. Let K 3d be some path connecting 
L 3d and L d inside the conformal quadrilateral R d d = Q d \ Q 3d (see Fig. 0). Then 

uj(v] K 3d ; Q d \ K 3d ) ^ const -u(v; L d \ Q d ) 

for some absolute positive constant. 

Proof. Note that u(v; L d ; Q d ) < u(v; L d ; Q d \ K 3d ) + u(v; K 3d ; Q d \ K 3d ). Thus, it is 
sufficient to prove that 

lu(v] L d ; Q d \ K 3d ) < const -u(v; K 3d ; Q d \ K 3d ). 

Furthermore, monotonicity arguments give 

u(v;L d ;tt d \K 3d ) «C / u(z; L d U L 3d ; Rf \ K 3d ) ■ uj(v; \dz\; tt 2d \ K 3d ) 

JL 2d 

and, in a similar manner, 

uo{v- K 3d - n d \ K 3d ) > / uo{z- K 3d - Rf \ K 3d ) ■ u(v; \dz\;Q 2d \ K 3d ). 
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Let L d = A d B d and so on (see Fig. HJ). Monotonicity arguments give 

u(z; K 3d ; R 3d \ K 3d ) > min {cu(z; A 3d A d ; R 3d ) , u{z; B d B 3d ; R 3d )} . 

Thus, it is sufficient to prove that 

u(z; A d B d U B 3d A 3d ; R 3d ) ^ const • min {cu(z; A 3d A d ; Rf) , u{z; B d B 3d ; R 3d )} 

for all z G L2 d . Due to the conformal invariance of harmonic measure, the last esti- 
mate follows from the uniform bounds on the extremal distances (conformal modulii 
of quadrilaterals) 

I 777, 

\A nd B nd B md A md (AndA md ; B md B nd ) ^ — log — , 1 n < m ^ 3. □ 

3.5. Convergence for the Poisson kernel normalized at a boundary point. 

Recall that HP 5 denotes the discretization of the upper half-plane EI = {z : Im z > 0} 
(i.e., the set of all faces, edges and vertices that intersect H, see Fig. EB)- As for 
bounded domains, denote by {a/ 5 }; Hp) the probability of the event that the 

random walk starting at u 5 G Hp first hits the boundary <9H r at a vertex x 5 G cfflf. 
It is easy to see (e.g., using the unboundedness of the free Green's function (12.31) or 
Proposition I2.1Q[ ) that 

£ w V;{^};Hr) = i- 

x s edw s r 

Let 

% 6 (u s ) = Imu s - Imx<5 ■oo 5 (u s ;{x 6 };C s + ) for u 5 G Hp . (3.8) 

x s £dW s r 

The function S* 5 is discrete harmonic in Hp, = on <9Hp and |3 <5 (m <5 ) — Im?/| ^.25 
for all u s G Hp (note that these conditions define uniquely). In particular, if 
Imu 5 G [35,5(5], then ^(u 5 ) x 5. Since is discrete harmonic and nonnegative, it 
implies 

Below we need some simple technical estimates for the discrete harmonic measure 
in the rectangle at a point lying near the boundary. Let o s G <9Hp be the boundary 
vertex closest to 0, 

R(S,T) = (S;S) x (0;T) C C, S,T > 0, 

Rf<(S, T) C T denote the discretization of this open rectangle, and 
L 5 r (S) = {v 5 G dR s r (S,T) : Inn/ 5 < 0} , 
U£(S,T) = {v 5 G dR s r (S,T) : Inn/ 5 ^ T} , 
V/(5,T) = {f 5 G dR^(S,T) : |Re?/| ^ 5} 

be the lower, upper and vertical parts of the boundary dR^(S,T) (see Fig. EB). 

Lemma 3.12. Let s ^ 2t > 0. If 5 is small enough, then 

u s (of nt ; U$(s, t); i? r (s, £)) x <f/t and c/(of nt ; V r 5 (s, i); i?f (s, £)) < const -5t/s 2 . 

Proof. We defer the proof until after that of Theorem 13.131 □ 
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From now on, we consider only discrete domains (fl*; a s ), a s G <9f2p such that 

R 5 r (S,T) cQ 5 r , L 5 r (S)cdQ 5 r and a 5 <£ L 5 r (S) (3.9) 

for some S, T > 0. All continuous domains (fl; a) appearing as Caratheodory limits of 
these (D, 6 ; a s ) satisfy 

R(S,T) C Q, [-IS; §S] C dQ and a £ \S}. (3.10) 

For continuous domains (fl;a) satisfying (13.101) . we define the Poisson kernel 
Po{-) = Po{-', a; Q) normalized at as the (unique) solution of the following problem: 

AP = in O, ^L Ua} = 0, Po^O and d y P (0) = 1. 
For discrete domains (f^a 5 ) satisfying (13.91) . we call = P s s (-;a 5 ;Q 5 r ) the discrete 



Poisson kernel normalized at o , if 

\dn s r \{a} 



A^=0 in ^, p*|_, vri =0 and P 5 AoL) = 5 (of, 



Note that 



PU-,a s ;Q s r )=u s (-,a s :Q s " ^ ( °* nt 



Theorem 3.13. T/ie discrete Poisson kernels P S s(-] a s ; f2f ) defined for the class of 
discrete domains satisfying (E2j) some 5*, T > are uniformly C 1 -close in Q s (see 
Definition \ 3. 6\) to the continuous Poisson kernels Po(-; a s ; fl s ), where Q s denotes the 
modified polygonal domain Q s with the part of the boundary L S (^S) C dQ s replaced by 
the segment [— ^S, ^S]. The rate of uniform convergence depends on S,T. 

Proof. The continuous Poisson kernel Pq is Caratheodory stable, so (13.51) holds true. 
Thus, it is sufficient to check (13.61) . Without loss of generality, we assume that S ^ 2T. 

Let (Q s ; u s ; a 5 ) ^{Q; u; a). Denote by T 5 G U$(±S, \T) the vertex closest to the 
point |T G C. Due to the discrete Harnack Lemma (Proposition 12.71 (ii)), the values 
of P£ s on UrioS, \T) are uniformly comparable with P^ S (T S ). Using Lemma [3.121 one 
obtains 

* ~ Pos(o 5 mt ) > const -P 5 oS {T 6 ) ■ o/(of nt ; U & r {\S, |D; RviiS, \T)) > const ~ • P & oS {T 5 ) 

and so the values P^ S (T S ) are uniformly bounded. Hence, applying the discrete Har- 
nack Lemma once more, one derives the uniform boundedness of P S s( 4 ', a 5 ; fif) on each 
compact subset of Q. Then, due to Proposition 13. 11 one can find a subsequence — > 
such that 

P s J(-,a Sk ;np)=tH 

uniformly on compact subsets of Q, where H ^ is some harmonic in Q function. We 
need to prove that H{u) = Po(u; a; Q) for each subsequential limit H. 

Repeating the arguments given in the proof of Theorem 13.101 one obtains that, 
first, for each r > the functions P s are uniformly bounded everywhere in flf n V s 
(in particular, everywhere in the rectangle R^(S,T)) and, second, H = on dQ \ {a}. 
Therefore, due to H ^ 0, 

H = fiPo(-; a; Q) for some /i ^ 0. 
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Let 

Q s (.) = P* s (.;a s ;n s r ) -&(■). 

By definition, the function Q s is discrete harmonic in R^(S,T), Q s = on the lower 
boundary L S T (S), Q s (of nt ) = and 

Q s (v) =4 ^Pq{v] a; Vt) — Imv as 5 = 5k — > 

uniformly on compact subsets of R(S,T). 

Since Po|[-s,s] = and (Po)' y (0) = 1, one has Po(x+iy) = y + 0(xy + y 2 ) as 
x, y — *■ 0. Thus, for any fixed s < S and t < T, 2t ^ s, the following hold true: 

Q 5 = on < const on ^(s,t), 

Q 5 =4 (/i-l)t + O(st) as 5 -> uniformly on E$(t, s). 

The normalization Q 5 (of nt ) = and Lemma [3.121 give 

\(ji-l)t + 0(st) + o(l) | •<?/< ^ const -5t/s 2 
as 5 —>■ 0. So, for any fixed s and t, one has \fi — 1| ^ const -(s + t/s 2 ). Thus, jM — l. □ 
Proof of Lemma 13.121 We consider two quadratic harmonic polynomials 

2\ „, „,2 Jl 



and 



(t+25) 2 \ y+25 x 2 -(y+2Sf 



h 2 (x+iy) = 1 + ^ 5-^- ■ ^— — + 



s 2 J t+25 s 2 

Their restrictions on F are discrete harmonic due to Lemma 12.21 (i). Moreover, it is 
easy to check that 

ht(v 5 ) < < h 2 {v s ), if v s G Lf(s) U V^?(s,t), 
< 1 < h 2 (v s ), if v 6 eU£(s,t). 

Hence, ^(V 5 ) c/fV 5 ; C/^(s, i); i?f (s, £)) < fc 2 (^) for all v 5 G R^{s,t). Thus, 

c/(u'; *);#f(M)) x 5/t, 
if, say, lmv s G [25,45] and | Reu 5 | ^ 45. Since u s is discrete harmonic, it implies 

co\of nt ;U^s,ty,R 5 r (s,t))-5/t. 

The upper bound for u s (v s ;V^(s,t); R^(s,t)) follows by the similar consideration of 
the quadratic harmonic polynomial 

= * + <„ + m + M-v) 

s 2 

which is nonnegative on Lf(s) U U[(s,t) and no less than 1 on V^(s, t). □ 
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A. Appendix 



A.l. Kenyon's asymptotics for the Green's function and the Cauchy kernel. 

Below we give a sketch of R. Kenyon's |Ken02] arguments with some small improvement 
concerning the estimates of the integrands for intermediate A. See also [Biick08j. 



Proof of the Theorem 12.51 Following J. Ferrand [Fer44j and Ch. Mercat [Mer07] 
we introduce discrete exponentials 



e(X } u;u ) = J[ - 



2k-l 1 . \, 

1 + o{U k - 



■U k -1 



j=0 - ^(Uk-Uk-x) 



(A.l) 



where P UQU = uoUiU2-..U2k-iU2k is a path from uq to U2k = u on the corresponding 
rhombic embedding (so u<ij G T and U2j-i G T* for all j). It's easy to see that this 
definition does not depend on the choice of the path. Since the angles of rhombi are 
bounded from and n, one can choose P UQU so that the following condition holds: 

for all j either (a) | arg(iij+i— Uj) — arg(w— u )\ < | or (b) Uj and Uj+2 are 



opposite vertices of some rhombus and | arg(«j +2 — Uj) — arg(w — u )\ < |. 
In particular, arg(u J+1 — Uj) — arg(n — u ) G (— for all j. 

Define (see R. Kenyon [Ken02j and A. Bobenko, Ch. Mercat and Yu. Suris [BMS05J) 

1 f logA 



G T (u;u ) 



]ir 2 i 



c 



A 



e(A, u; uq) dX. 



(A.2) 



where C is a curve which runs counter clockwise around the disc of (large) radius R from 
the angle arg (u— u ) — n to arg (u— uq) + tc, then along the segment e ta,TS ^ u ~ u °^[—R, — r], 
then clockwise around the disc of (small) radius r and then back along the same segment 
(the integral does not depend on the log branch, since e(0, u; u ) = e(oo, u; u ) = 1). 

This function is discrete harmonic away from Uq since all discrete exponentials 
are harmonic (as functions of u) and one can use the same contour of integration for 
all u s ~ u. Furthermore, G(u ; u ) = and, by straightforward computation, 

Gr(u s ;u ) = 6>sCQt ^ ; jf -u s ~ it, so fx s r (u ) ■ A s G r (u ;u ) 



1. 



7T 



Rotating and scaling the plane, one may assume that arg(w — uq) = and 5=1. 
It's easy to see that the contribution of intermediate A = — t < to the integral in 
(IA.2[) is exponentially small. Indeed, in case (a) one has 

2 



1 



2 C 



l + -e 



1 



8t cos Pj 



t 2 + At cos (3j +4 



exp 



8t cos Pj 



where Pj 



a.rg(uj + i — Uj ) and cos Pj = 
[1 - |e^)(l - |e ift+i ; 



Re(uj + \ — Uj) > 0. Similarly, in case (b), 



l + |e^)(l + |e 



'i+i 



exp 



8t(cos Pj + cos Pj+i] 



H + 2Y 



due to cos/3j+cos/5 J+ i = Ke(uj + 2—Uj 



|e(— t, u; u ) I ^ exp 



> 0. Thus, 

At(u- 



■uo) 



(t + 2) 
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and the asymptotics of flA.2[) as \u — u \ — > oo are determined by the asymptotics 
of e(X,u;uo) near and oo. Some version of the Laplace method (see [Ken02j and 
|Biick 08j) gives 

7Euler + log 2 



G r (u;u ) 



2tt 



log \U — Uq\ + 



2tt 



+ O(\u—u \ 



where the remainder is of order \u— Uq\ 2 due to 



dX 2 



log(e(X, u;u )) 



d 2 



log(e(A,«;ito)) 



0. 



A=oo 



a=o dX 2 

The uniqueness of G(-;uq) (and G = ReG) easily follows by the Harnack inequality 
(Corollary 12. 81) . Indeed, the difference G = Gi(-;uq) — G2(-;uq) is discrete harmonic 
on T and max| u j^ R \H(u)\/R — > as u — > oo, so G(u) = G(uq) = for all u. □ 

Proof of Theorem 12.141 As in Theorem l2.5l K(-; zq) can be constructed using (mod- 
ified) discrete exponentials. Similarly to flA.ll) denote 

1 



e(X,v j ;z ) 



J 



1,2,3,4, 



for the vertices of the rhombus v^v^v^v^ centered at z and, by induction, 



(A.3) 



e(X,w;z ) 1 + %(w-u) 

- f- for all w ~ u, u G T, wGr 

1 - ~(iu-u) 



e(A,u;2 ) 

As above, all e(A, ■; zq) are discrete holomorphic on A. Define (see [Ken 02j) 



K(v; Zq) = - 

71 



(v— «o)oo 



e(A, v; z )dX, 



where the integral being, say, along the ray arg£ = arg(t> — zq) ± n (taking the path 
from zq to u as in the proof of Theorem 12.51 one guarantees that all poles of e(-, v; Zq) 
are in | argA — arg(t> — zq)\ < it). Then K(-; zq) is holomorphic everywhere except zq. 
The straightforward calculations give 



Vzov? ■ K{v°a z Q ) 



arg(uj -u}^) 



arg(vJ-i;J +1 ) 



so 



7T 



4(^o)-[^ir](z ;^ 



1. 



Let 5=1. As in Theorem 12. 5[ the integrand is exponentially small for intermediate A. 
One has 

e(X,v; zq) = exp [X(v-z ) + O(\X\ 2 \v-z \)] , A -> 0, 
and 

'|u-Zo| 



4r 2 

e(A,v; zb) = — • exp 



A 



+ 



IAI 



A 



oo, 



/ (!) ( 2 )\ • (■ C 1 ) ( 2 h _n 

where r = e jargin ° j , if v G T, and r = e iarg[ - w ° w o >, if v G T* (the factor comes 



from (1A.3j) ). Summarizing, one arrives at 



K(u; z ) 



71 



1 T 2 
+ 



V-Zq V-Zq 



+ 



V-Zq\ 



Pr 



7T 



-; r 



+ 



Finally, K(-,zq) is unique due to Corollary 12.81 



□ 
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A. 2. Proof of the discrete Harnack Lemma. Below we recall the modification 
of R. Duffin's arguments [Duf53| given by U. Bucking |Biick08]. For the next it is 
important that the remainder in ( 12.31) is of order 5 2 \u— uq\~ 2 (and not just 5\u— mo| _1 ). 

Proposition A.l. Let Uo&T and 5. Then 

u 5 (u ; {a}; B s T (u Q , R)) x - for all a £ dB s r (u , R), 

R 

i.e., consti -5/R ^ oj 5 {uq] {a}; fip) ^ const2 -S/R for some positive absolute constants. 
Proof. One has R ^ \a — uq| ^ R+25 for all a E dB^(uo, R). Therefore, (12.31) gives 



n ( n 1 i \ u ~ u v\ 



5 ( 5 2 5 2 \ 
^ — ^ + const ■ i — + — 



nR \\u—uo\ 2 R 2 J 

for all u «o- In particular, if R/S is large enough, then 

8 



\G B 6( UQjI q(u;Uo) \ x — for all u G B T (u , R) : R—55 ^ \u — « | ^ R—35. 

Since Green's function is discrete harmonic and nonpositive near dB$,(u ,R), the same 
holds true for all a int : (a; Oi nt ) e dB^(u ,R). In view of (I2.5p . this gives the result 
for sufficiently large R/5. For small (i.e. comparable to 5) radii R the claim is trivial, 
since the random walk can reach a starting from Uq in a finite number of steps. □ 

Proposition A. 2 (mean value property). Let H : By(uo,R) — ► M be nonnegative 
discrete harmonic function. Then 



7TR 2 

., I. .4- 



jieIntB^(Mo,-R) 



. , SHjup) 

^ const — . 

R 



Proof. Let 



F[u) = G r {u; m ) ^j - + ' u r = 5 rK, 

Note that [A 5 F](w) = -(vri? 2 )- 1 for all u ^ m (see LemmaOJi)). Using (Q, it is 
easy to see that F{u) = 0(S 2 /R 2 ), if \u — R\ ^ const -5. The discrete Green's formula 
gives 

%)-^ H i u )^{u) = J2 ^n6 aaint -[F ± (a)H(a iQt )~F ± (a int )H(a)], 



uSlnt B r 



where the functions F^ = F ± const -5 2 / R 2 are positive/negative, respectively, near 
dBf(u , R). Using H ^ 0, one obtains 

5 2 1 5 2 

-const— 22 H ( a ) ^ H ( u o) ~ ^5 /J ^(^)«r(«) ^ const ■— H(a int ), 

Both sums are comparable to 5H{uq)/R due to Proposition lA.il □ 
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Proof of Proposition YI7f\ (i) Applying the mean value property for B^(uq,R) and 
B^(ui, R—28) and taking into account H{u\) ^ const •if(wo), one obtains 

7rR 2 H(u ) -tv(R-25) 2 H( Ui ) = H(u)fi s r (u) + O(5RH(u )). 

u£Bf,(u ,R)\B*,(u u R-25) 

Proposition IA.1I gives 

H(u)jj s r (u) x V 8 2 H(u) = O(8RH(u )), 



E 



«£Bf (u ,R)\Bt, (m ,R-25) ueB* (u ,i?)\B£ (m ,R-2S) 

so R 2 (H( Ul )-H(u )) = 0(6RH(uq)). 

(ii) Let Ui = v viV2...Vk-iVk = u 2 be some path connecting u\ and u 2 inside Bp(u ,r) 
(one can choose a path such that k ^ const -r/8). Since B^ Vj ,R-r) C B^{u ,R), 

fe-i 



#(u 2 ) _ ^ 



n 



1 + const ■ 



5 



R-r 



- const -r/8 


r 


^ exp 


const ■ — 

R — r 



□ 
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